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Abstract 


A  one-dimensional  fluid  model  for  a  coupled  system  of  klystron  output 
cavities  is  developed.  The  model  is  exercised  for  both  realistic  and  opti¬ 
mized  beam  current  modulation  waveforms,  and  the  rf  conversion  effi¬ 
ciency  is  calculated.  An  optimized  second  harmonic  velocity  modulation 
and  associated  current  modulation  are  calculated  and  inputted  into  the  cav¬ 
ity  simulation.  Optimum  bunching  lengths  and  rf  conversion  efficiencies 
are  determined. 
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1.  Introduction 


Klystron-like  (klystron,  multibeam  klystron,  reltron,  RKA)  rf  (radio  fre¬ 
quency)  sources  are  very  suitable  for  high-peak  and  average-power 
microwave  applications,  including  directed  energy  warfare  (DEW).  This 
suitability  is  due  to  a  configuration  in  which  the  electron  beam  collector  is 
separate  from  the  rf  generation  cavities  (unlike  the  magnetron  or  crossed 
field  amplifier  (CFA)).  Because  of  this  separation,  heating  in  the  rf  cavities 
is  due  only  to  conduction  currents  flowing  in  the  cavity  wall  and  not  to 
electrons  impinging  on  the  walls  of  the  rf  cavities  (if  the  tube  is  correctly 
designed).  The  major  deficiencies  of  klystron-like  rf  sources  (compared  to 
magnetrons  and  CFAs)  are  a  relatively  high  operating  voltage  and  a  sig¬ 
nificantly  lower  efficiency.  In  this  report,  I  address  the  problem  of  high 
voltage  in  a  klystron-like  device  by  using  a  high  perveance  (>10  /xperv) 
electron  beam  as  is  done  in  the  multibeam  klystron  amplifier.  The  prob¬ 
lem  of  increasing  the  klystron  efficiency  centers  on  the  particular  current 
modulation  impressed  upon  the  electron  beam.  To  determine  an  achievable 
electron  beam  modulation  that  will  maximize  the  klystron  efficiency  while 
minimizing  the  number  of  cavities  and  hence  the  length  of  the  klystron,  I 
have  developed  a  fast  one-dimensional  fluid  simulation  of  a  klystron  out¬ 
put  cavity.  This  simulation  will  allow  me  to  perform  a  comprehensive  para¬ 
metric  study  on  the  beam  modulation  required  to  maximize  the  rf  conver¬ 
sion  efficiency  of  the  klystron.  After  the  optimum  current  modulation  is 
determined,  a  realizable  velocity  modulation  that  will  produce  an  approx¬ 
imation  to  the  optimum  current  modulation  is  derived,  and  the  rf  conver¬ 
sion  efficiency  for  such  a  current  modulation  is  determined  with  the  fluid 
simulation.  In  particular,  the  optimum  ballistic  bunching  distance  for  the 
derived  velocity  modulation  is  determined  with  the  fluid  simulation  from 
a  parametric  study. 
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2.  Klystron  Output  Cavity  System 


2.1  Multiplet  Output  Cavity  System 

The  purpose  of  the  klystron  output  (or  multiplet  output)  cavity  system  is 
to  efficiently  convert  the  kinetic  energy  in  the  current-modulated  electron 
beam  into  rf  power.  An  example  of  an  output  cavity  system  that  is  analo¬ 
gous  to  the  doublet  buncher  cavity  [1]  is  shown  in  figure  1.  The  only  dif¬ 
ferences  in  the  configuration  of  the  buncher  and  the  output  cavities  is  that 
in  the  output  cavity,  no  generator  is  in  series  with  the  load  impedance  and 
the  electron  beam  entering  the  cavity  is  highly  current  modulated. 

As  in  a  previous  paper  that  analyzed  the  buncher  cavity  [1],  I  will  analyze 
the  general  case  of  a  multiplet  output  cavity  system.  Unlike  the  buncher 
cavity,  a  single-harmonic  time  analysis  would  not  be  useful  for  the  output 
cavity.  First,  because  of  the  nature  of  ballistic  focusing,  the  current  modu¬ 
lation  of  the  electron  beam  entering  the  output  cavity  system  is  very  rich 
in  higher  harmonics.  Second,  since  current  gains  (over  the  current  injected 

Figure  1.  Doublet  . . .  _ 

output  cavity.  (Variable 

will  be  explained  later  _  L 

in  report.) 
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into  the  buncher  cavity)  of  greater  than  20  are  achievable,  the  peak  current 
in  the  output  cavity  system  could  approach  the  static-limiting  current  for 
the  cavity  gap.  In  the  output  cavity  system,  the  effects  of  both  space-charge 
and  higher-order  cavity  modes  must  be  considered.  The  space-charge  ef¬ 
fects  may  increase  the  efficiency  of  the  output  cavity  by  transferring  energy 
from  the  higher-current  harmonics  to  the  fundamental  frequency. 


2.2  Driven-Mode  Equations 

In  the  output  cavity  system,  the  current  distribution  driving  the  cavity  is 
not  guaranteed  to  be  sinusoidal  in  time.  To  account  for  this  situation,  I 
must  use  the  general  expression  for  the  amplitude  c^k(t)  of  the  Ith  driven 
mode  in  the  kth  cavity  of  an  N  cavity  output  cavity  system.  The  following 
equation  (1)  is  the  same  as  equation  (7)  in  a  previous  paper  [1],  but  with 
Vgen  =  0,  Zgeu  replaced  by  Z\oa(\,  and  the  gap  and  loop  voltages  in  each  cav¬ 
ity  being  the  sum  of  the  gap  and  loop  voltages  of  the  fundamental  cavity 
mode  (l  =  0)  and  the  higher  harmonic  modes  (1  <  /  <  Nk)  in  the  kth  cavity: 


$,*(*)  +  ulkchk(l)  =  — t — — - < 

e0  dV  |4fe(xfe)|2 

Now  I  introduce  the  phase  variable  ip  =  u>t,  where  u  is  the  drive  frequency 
of  the  buncher  cavity.  While  I  cannot  expect  c^k{t)  to  be  sinusoidal  in  time, 
I  can  assume  that  it  will  achieve  a  periodic  equilibrium  with  period  2-k/uj. 
With  this  change  of  variable,  let  zk  =  zk/Dk  and  u jkk  =  u>ik/ u,  so  that  the 
driven  mode  equations  become 
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The  system  of  J^k=iNk  equations  described  by  equation  (4)  is  of  a  second 
order  and  thus  is  not  suitable  for  a  solution  (when  coupled  to  the  electron 
fluid)  by  standard  ordinary  differential  equation  solvers  (Runge-Kutta,  etc). 
I  will  convert  equation  (4)  to  a  system  of  2 J2k=i^k  first  order  differential 
equations  by  introducing  the  variables  Vu--  These  variables  are  the  driven 
gap  voltages  for  the  /th  mode  of  the  fcth  gap  in  the  multiplet  output  cavity 
system.  The  relationship  between  Vik  and  c/,*  is  given  by 

YlA  - _ l  Yjjj  gaP  A.  .  <*\  ic\ 


VIA  Dk- 


- Cl.k(t )  • 


The  additional  J2k=i  ^a-  first  order  equations  required  to  complete  the  sys¬ 
tem  are 

-JZ  (T/.A-gapQ.A-)  =  IVlAVlA  •  (6) 

Using  equation  (6)  to  eliminate  the  second  derivative  from  equation  (4) 
gives 
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Now  I  introduce  the  dimensionless  variables  V/j .  and  ckk  defined  by 


Via  = 


.  Via  gap 

Q.A-  =  t—  Q.A- 
Miorm^ 
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so  that  the  dimensionless  first  order  driven-mode  equations  become 
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dip 
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(12) 


I  shall  leave  the  normalization  constant  ynormfc  undefined  for  now  so  that  I 
can  select  a  value  for  it  later  that  will  lead  to  the  greatest  simplification  of 
the  fluid  equations.  Note  that  for  a  lossless  cavity,  <1^,  A and  G/  j.  will 
be  real  and  the  system  of  equations  describing  the  electron  beam  cavity 
system  will  be  real. 


2.3  Space-Charge  Fields 


For  the  Zcth  planar  gap,  Poisson's  equation  is 

d  (Pk _ _ pk  jzki  p) 

dzl  ~  e0 

The  same  equation  with  a  normalized  position  variable  (%)  becomes 

<Ppk  =  DlpkjzkPP) 

dzl  £o 


(13) 


(14) 


with  boundary  conditions  <pfc(0)  =  <p\.  (1)  =  0.  Equation  (14)  can  be  directly 
integrated,  giving 
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Equation  (20)  gives  the  contribution  of  the  space  charge  electric  field  to  the 
total  electric  field. 


Fluid  Equations 

The  fluid  equations  for  the  A  th  planar  cavity  gap  (/th  cavity  mode  excited) 
are 


dpk 

dt 


+  pk 


du^ 

dzk 


+  u'~k 


dpk 

dzk 
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and 


+  u- 
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Transforming  equations  (21)  and  (22)  into  ones  that  will  use  the  dimension¬ 
less  variables  p  and  zk  gives 
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Substituting  equation  (20)  into  equation  (24)  gives 
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Making  equation  (23)  dimensionless  requires  that 


(25) 
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and  then  equation  (23) 


ucLumo 
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dp  dzk  dzk 


(27) 


The  form  of  equation  (27)  allows  pk  to  be  normalized  in  any  way  to  simplify 
equation  (25).  First  I  write  equation  (25)  with  the  chosen  u~k  normalization 
(eq  (26))  to  show  how  to  normalize  pk: 
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By  inspecting  equation  (28),  I  see  the  correct  normalizations  are 

W  =  and 
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These  equations  transform  equations  (27)  and  (28)  into 
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A  normalized  load  impedance  Zik  is  defined  by 
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so  that  equation  (12)  becomes 
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and  the  system  that  must  be  solved  for  determining  the  performance  of  the 
output  cavity  consists  of  equations  (11),  (31),  (32),  and  (34). 

2.5  Finite  Element  Expansion 

To  solve  the  output  cavity  equations  with  the  use  of  finite  element  method, 
the  first  step  is  to  expand  uZk(zk,  p)  and  &k{zk-  p)  as  a  set  of  at  least  once 
differentiable  axial  basis  functions  ( £/, ) .  Note  that  in  subsequent  formu¬ 
las,  NZk  denotes  the  number  of  axial  basis  functions  for  the  kth  gap.  The 
basis  functions  must  be  differentiable  at  least  once,  since  I  am  solving  a 
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dVi,k 

dip 


system  of  first-order  equations.  I  am  also  restricting  the  basis  functions  to 
be  real  to  speed  up  the  numerics  of  the  solution  process.  The  finite  element 
expansions  for  the  velocity  and  charge  density  fields  are 

Uzt(Zk,V>)  =  5Z  “**•».  Sm  (**)  and  (35) 

m  —  1 
N-k 

Pk(Zk,<P)  =  ^Pk.m  (?)  (4)  ■  (36) 

771  =  1 


In  the  following  equations,  prime  when  applied  to  a  function  will  denote 
differentiation.  Second,  I  define  the  auxiliary  basis  functions  as 


tm  (zk)  —  [  dz[.  S*  (4)  - 
Jo 

and  the  dyadic  matrixes  as 
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(40) 

(41) 
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Then  the  third  step  is  to  reduce  equation  (34)  by  noting  that  the  finite  ele¬ 
ment  expansion  allows  me  to  substitute 
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giving 


A’  A'n 


djl.kh.k  ~  &d>l.k  X  (  G/.*  £  Ef41 + I  •  <«« 

77?  =  1  71=0  7?*?7?  771  =  1  71  =  1 


This  is  the  only  effect  of  the  finite  element  expansion  upon  the  cavity  equa¬ 
tions.  Now  I  will  reduce  the  charge  conservation  equation  to  its  finite  ele¬ 
ment  equivalent.  First  I  rewrite  equation  (31)  as 


dpk  djihjk) 

dp  dzk 


(45) 
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Next  I  multiply  equation  (45)  by  sf  (zk)  and  integrate  from  0  to  1  by 

<46> 


Substituting  equations  (35)  and  (36)  and  integrating  by  parts  give 
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Using  the  definitions  for  and  Cfmn  reduces  the  finite  element  approxi¬ 
mation  of  the  charge  conservation  equation  to 


E  AtJ^  =  E  £cL,k,m^.,  +  4'  (0)  Jk  (0,  v)  =  0  ■  (48) 

m=l  ^  m=l  n=l 

The  reduction  of  equation  (32)  (fluid  equation)  proceeds  in  a  like  manner: 


14, fc 


(49) 


Treating  each  part  of  equation  (49)  separately  gives  the  following  finite  ele¬ 
ment  expansions: 


JV, 


Y'  Ak 

^  1171  dip 


m=l 


N*k  N*k 


~  Y  YClmn^k,A,n  ~  ~sf  (0)  UZk  (0,  pf  , 

771=1  77  —  1 

/  Nk  ~ 

=  B\  5>n,* 

\  n=0 


[  dzksi  [  dz  f  dz"  Pk 
Jo  Jo  Jo 


Y  DlmPk,m  ,  and 

m=  1 
"•k 

Bl  Ys  ^rnPk.m  ■ 

777  =  1 


(50) 

(51) 

(52) 

(53) 


Substituting  equations  (50),  (51),  (52),  and  (53)  into  equation  (49)  gives  the 
finite  element  form  of  the  fluid  equation. 
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k  dii 

ST  a  k  auzkj 

lm  dip 
m~  1  ^ 


-i  JV.fc 

=  nE  ECton%. 


m=l  r?~  1 


+£«*  (°)  «** 


EAmPfc.m-  (54) 

777—  1 


Finally,  the  boundary  conditions  on  the  beam  velocity  and  current  density 
fields  are  given  by  equations  (55)  and  (56)  for  k  >  1: 

AVi 

“**(' 0,(p)  =  £  (l)uSk  k_,  m(p)  and  (55) 

777  —  1 

nh-, 

Jk  (0,  p)  =  UZk  (0.  p)  Y,  (1)  Pk.k~l.ni  (p)  -  (56) 

777  —  1 


For  A:  =  0,  the  beam  modulation  functions  uZo  (0,  <p)  and  J0  (0.  p)  are  either 
assumed  or  calculated  from  the  ballistic  bunching  equations. 


2.6  Basis  Selection  for  Finite  Element  Formulation 


The  basis  functions  selected  for  the  finite  element  formulation  should  allow 
for  efficient  evaluation  of  the  derivatives  in  equations  (11),  (44),  (48),  and 
(54)  so  that  a  standard  differential  equation  solver,  such  as  Runge-Kutta, 
can  be  used  to  solve  the  system.  Additionally,  the  basis  functions  should 
be  in  a  form  such  that  relatively  few  functions  are  required  to  approximate 
the  velocity  and  charge  distributions.  Since  some  of  the  current  modula¬ 
tion  waveforms  that  will  be  used  to  model  the  electron  beam  injected  into 
the  output  gap  have  large  derivatives,  the  basis  functions  must  be  able  to 
model  charge  densities  that  have  a  rapid  spatial  variation.  Because  of  this 
property,  basis  functions  consisting  of  e±27n/5  were  not  considered,  since 
large  values  of  /  would  be  required  to  model  rapidly  varying  functions.  In 
the  end,  I  selected  the  basis  set  for  the  piecewise  linear  functions.  If  the  in¬ 
terval  [0. 1]  of  the  A-th  output  gap  is  partitioned  into  N~k  intervals  by  the 
grid  {z,k0, ...  ,zkt. ... , ZkNt  }/  where  zk„  —  0  and  zkN_  =  1,  the  NZk  piece- 
wise  linear  basis  functions  defined  by  this  grid  are 


si 


-  zA-  o  <  4  <  zk  j 
zh 

0  <  **  < 


(57) 
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Si  (. 4 ) 


sNZk  (4) 


0 

4  ~  4t-i 
4,  -  Zklzl 
hi+i  ~  4 

4(+1  -  4, 

0 


0  <  4  <  4,_i 
4,_!  <  4  <  4, 

h,  <  4  <  ht+1 
hl+1  <  h  <  l  , 


and 


0 

1-4 
1  - 


0  < 

i 

► 

4JV^_1  <  4  <  l 


(58) 


(59) 


so  that  an  arbitrary  function  fk  (4)  that  takes  on  the  values  fk  (zkl)  =  ff 
in  the  /;:th  output  gap  may  be  represented  by 


/fe(4)  =  E/^(4)  •  (60) 

i=i 

This  choice  of  basis  results  in  sparse  arrays  for  and  Cfmn  and  a  trian¬ 
gular  array  for  D^n.  All  arrays  are  evaluated  using  a  C++  class  library  for 
piecewise  polynomial  functions.  First  the  basis  set 

{sf(4),  1<1<JV**}  (61) 

is  defined  as  an  array  of  piecewise  linear  functions.  This  array  of  functions 
is  then  integrated  (one  of  the  standard  operations  for  the  class)  to  generate 
the  associated  array 

{t?  (4) ,  1  <l<NZk}  (62) 

of  piecewise  quadratic  functions.  Then  the  functions  in  the  arrays  are 
appropriately  multiplied  to  form  the  piecewise  polynomial  functions 
(quadratic  and  cubic)  that  are  the  integrands  in  equations  38  through  42. 
Multiplication  of  piecewise  polynomial  functions  is  also  a  standard  class 
operation.  The  integrals  are  then  analytically  evaluated  (another  standard 
class  operation).  The  piecewise  polynomial  class  library  simplifies  evalu¬ 
ating  the  required  arrays  for  any  set  of  basis  functions  that  are  piecewise 
polynomials. 

From  a  computational  point  of  view,  the  selected  basis  set  is  very  attrac¬ 
tive  since  the  array  Aknn  is  a  tridiagonal  matrix.  This  property  allows  rapid 
evaluation  of  the  derivatives  in  the  finite  element  equations  with  the  use 
of  standard  algorithms  [2,  p  50],  even  though  the  basis  functions  are  not 
orthogonal  (Ajj,n  is  not  diagonal). 
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Once  I  selected  the  basis  functions  and  evaluated  the  finite  element  ar¬ 
rays,  the  next  step  was  to  solve  the  system  of  coupled  nonlinear  differential 
equations  defined  by  the  equations  (11),  (44),  (48),  and  (54).  The  method  se¬ 
lected  was  the  Cash-Karp  variation  of  the  Runge-Kutta  method,  with  adap¬ 
tive  step  size  control  [2,  p  717].  This  method  is  very  robust,  can  deal  with 
rapidly  varying  functions,  and  will  optimize  the  step  size  of  the  integra¬ 
tion.  This  method  was  implemented  as  a  C++  virtual  base  class,  with  the 
derived  class  supplying  the  required  derivative  functions  to  the  differential 
equation  solver  base  class. 


2.7  Calculation  of  rf  and  Beam  Kinetic  Power 


The  whole  purpose  of  the  klystron  is  to  produce  rf  power.  The  output  cavity 
is  where  the  rf  is  extracted  from  the  klystron.  To  calculate  the  rf  power 
extracted,  one  needs  to  calculate  the  load  current,  which  is  also  the  loop 
current,  I\oop  (t).  This  is  given  by 


-floop  (f)  = 
-fioop  (y)  = 

-fioop  (y)  = 

The  rf  power  averaged 


/  j  /  J  t  *  n.k  loop  • 


^load 


A=1  n—0  Un-k 
N  Nk 


Y'  1  dc^ 

Z^  Z^,~,  .  A  „  vn.k  loop 


^load  n^o^n.k  dp 

r2  N  n.  \  2  Nk  f. 


me*  y-  fu)Dky  y  VnM 

eZ^h  V  C  /  n^0Gl,  k 

over  one  rf  cycle  is  then 


and 


(63) 

(64) 

(65) 


(Ph)  =  fUdt  /i„op(t)2  and  (66) 

Jo 

(Prf)  =  ^ff^dpIloop(p)2  •  (67) 

Since  the  load  current  can  have  significant  harmonic  content,  I  performed 
a  Fourier  expansion  of  7ioop  (p)  and  calculated  the  power  associated  with 
each  harmonic  to  determine  the  harmonic  power  spectrum  of  the  klystron 
for  each  drive  frequency,  u. 


Since  I  used  approximate  models  to  determine  the  behavior  of  the  out¬ 
put  cavity  system,  I  needed  to  implement  diagnostic  tests  to  determine  if 
the  fundamental  conservation  laws  are  being  obeyed  by  our  approximate 
model.  To  do  this,  I  needed  to  calculate  the  average  electron  beam  kinetic 
power  at  the  exit  of  the  output  gap.  The  kinetic  power  density  of  a  one¬ 
dimensional  fluid  is  given  by 


dP\,(.  1  ( m\  3 

~dS~  ~  2  \  e)  PzU'~  ' 


(68) 
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For  a  beam  of  radius  Rb,  the  instantaneous  beam  kinetic  power  at  the  out¬ 
put  of  the  fcth  output  gap  is 


^ke(l)  =  (7)  =  Wo  {k°Rb) 2  {k°Dk) 3  (^r)  P^zk(l)3  >  (69) 


where  /cq  =  u>/c.  Define 


-  0.69313  GW 


(70) 


as  the  power  normalization.  The  average  beam  kinetic  power  during  one  rf 
cycle  is  then 

2D  /*27T 

(PL(l))  =  -f(k0Rb)2(koDkf  /  dypk{hv)uZk{l,V?  ■  (71) 

4  Jo 


With  the  use  of  the  normalizing  power  defined  in  equation  (70),  the  rf 
power  becomes 


PZTT  {  1 '  K 

/  d?  E(W2  E; 

JO  l  7)  n=0 


and  the  conservation  law  is 

(Fk°e(0))  -  Jim  {(Pa)  +  (P£(l))}-  (73) 

Equation  (73)  states  that  after  oscillatory  equilibrium  is  reached,  the  aver¬ 
age  kinetic  power  of  the  beam  entering  the  multiplet  output  cavity  is  equal 
to  the  sum  of  average  rf  power  to  the  multiplet  output  cavity  load  and 
the  average  kinetic  power  of  the  beam  exiting  the  multiplet  output  cavity. 
This  is  a  statement  of  conservation  of  energy.  Averages  over  equilibrium  rf 
cycles  must  be  used  because  of  the  energy  stored  in  the  cavity  and  space 
charge  fields. 


2.8  Output  Cavity  Operation  with  Specified  Current  Modulation 

After  being  coded,  the  fluid  simulation  of  the  output  cavity  was  run  for 
the  modulation  and  electron  beam  parameters  shown  in  table  1.  The  bal¬ 
listic  current  modulation  referenced  in  table  1  is  for  a  standard  two-cavity 
klystron,  where  a  buncher  gap  imparts  a  small  (5  percent)  velocity  mod¬ 
ulation  to  the  electron  beam,  and  the  beam  is  then  drifted  over  a  distance 
that  will  provide  optimum  beam  bunching.  The  equations  describing  this 
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Table  1.  Summary  of 
output  cavity  simulation 
results  for  10-pperv 
electron  beam.  No  ballistic 
drift  length  associated 
with  25  percent  gate 
modulation.  A  5  percent 
velocity  modulation  is 
assumed  for  a  ballistically 
generated  current 
modulation. 


Modulation 

Drift  length 
(cm) 

Beam 

power 

Voltage 

(kV) 

Current 

(A) 

Efficiency 

Ballistic 

25 

20  kW 

5.25 

3.8 

0.560 

Ballistic 

63 

2  MW 

33.14 

60.3 

0.565 

25%  gate 

— 

20  kW 

5.25 

3.8 

0.894 

25%  gate 

— 

2  MW 

33.14 

60.3 

0.896 

type  of  current  modulation  and  the  numerical  equivalent  of  the  treatment 
[3,  p  300]  are  given  in  section  3.1  of  this  report.  The  25  percent  gate  current 
modulation  is  simply  imposed  on  the  electron  beam  entering  the  output 
cavity  with  the  peak  beam  current  normalized  so  that  the  average  beam 
power  for  both  types  of  current  modulation  is  the  same.  The  modulation 
is  trapezoidal  in  time  with  rise  and  fall  times  that  are  5  percent  of  the  full- 
width-half-maximum  (FWHM)  duration  of  the  pulse  while  the  FWHM  du¬ 
ration  is  25  percent  of  an  rf  cycle.  For  all  the  cases  treated  in  this  report,  the 
rf  frequency  is  1  GHz  with  an  rf  cycle  time  of  1  ns.  Once  I  selected  the  mod¬ 
ulation  and  the  beam  parameters  (voltage  and  perveance),  I  ran  the  code 
for  several  instances  of  beam  loading  (combination  of  gap  gain  and  load 
impedance).  For  each  run,  the  system  was  stepped  forward  in  time  until 
it  achieved  oscillatory  equilibrium  or  until  the  electron  velocity  field  be¬ 
came  negative  at  some  location  leading  to  electron  trajectory  crossing.  This 
criterion  was  used  because  the  fluid  model  was  only  applicable  if  the  elec¬ 
tron  trajectories  within  the  buncher  cavity  gap  did  not  cross.  If  oscillatory 
equilibrium  was  achieved,  the  beam  loading  was  increased  until  oscillatory 
equilibrium  was  no  longer  possible;  then  the  beam  loading  was  slightly  de¬ 
creased  until  oscillatory  equilibrium  returned.  This  value  of  beam  loading 
is  considered  to  be  optimum — it  results  in  the  maximum  efficiency  for  the 
conversion  of  beam  kinetic  energy  to  rf  power. 

A  summary  of  the  code  runs  is  shown  in  table  1  with  the  detailed  results 
displayed  in  figures  2  through  9.  Figures  2, 4,  6,  and  8  show  both  the  input 
beam  current  and  velocity  modulations  ( It,  (0)  and  0Z  (0)  in  the  small  graph 
on  the  right  side  of  the  figures),  the  average  (over  one  rf  cycle)  rf  power 
(Pr f)  and  beam  kinetic  power  and  the  minimum  value  of  the  velocity 
field  in  the  cavity  gap  during  the  rf  cycle,  min  (/?-).  Currently,  the  cavities 
are  assumed  to  have  a  simple  transmission  line  geometry  as  in  figure  1  so 
that  the  eigen  modes  can  be  calculated  with  a  simple  analytic  expression. 
Figures  3,  5,  7,  and  9  show  the  detailed  behavior  of  the  beam  current  P  (1) 
and  beam  velocity  /3Z  (1)  modulation  at  the  output  of  the  cavity  gap.  Ad¬ 
ditionally,  the  voltages  across  both  the  load  Vjoac)  and  cavity  gap  Vgap  are 
shown. 
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The  most  striking  of  the  results  shown  in  table  1  is  that  a  simple  25  percent 
duty  factor  trapezoidal  beam  current  modulation  produces  a  significantly 
higher  conversion  of  electron  beam  kinetic  power  to  rf  power  than  the  stan¬ 
dard  ballistic  electron  beam  current  modulation.  This  effect  is  independent 
of  beam  power.  The  2-MW  beams  have  the  same  conversion  efficiencies  as 
the  20-kW  beams  with  the  same  modulation  waveforms.  The  explanation 
of  this  effect  can  be  found  by  comparing  the  beam  current  and  velocity 
fields  in  figures  3  and  7.  In  figure  3  when  (n/2  <  <p  <  37r/2),  the  current  at 
the  cavity  output  is  small,  and  the  beam  velocity  is  large.  When  the  velocity 
is  small,  the  beam  current  is  large.  However,  if  I  compare  the  beam  current 
under  similar  conditions  as  shown  in  figures  3  and  7,  the  beam  current  in 
figure  7  is  much  closer  to  zero  for  large  values  of  beam  velocity  compared  to 
the  beam  current  in  figure  3.  The  modulation  used  in  figure  7  gives  a  much 
higher  efficiency  than  the  standard  optimum  ballistic  modulation,  because 
the  modulated  beam  current  entering  the  output  cavity  is  zero  for  almost 
75  percent  of  an  rf  cycle.  For  the  trapezoidal  modulation,  all  the  electrons 
in  the  beam  slow  down  and  give  up  energy  to  the  rf  field.  For  the  ballistic 
modulation,  the  beam  current  is  never  truly  zero  and  a  significant  number 
of  the  electrons  entering  the  output  cavity  will  have  a  phase  that  causes 
them  to  be  accelerated  by  the  rf  field  in  the  cavity,  reducing  the  rf  output 
power.  To  increase  the  efficiency  of  the  klystron,  I  must  decrease  the  modu¬ 
lated  beam  current  for  those  phases  that  allow  acceleration  of  the  electrons. 
The  next  section  of  this  report  develops  a  method  for  implementing  such  a 
modulation. 

Before  moving  on,  one  numerical  aspect  of  the  fluid  simulation  requires 
some  commentary.  Note  that  on  the  trailing  edges  of  the  modulated  current 
/{,  (0)  exiting  the  output  cavity  (figs.  3, 5, 7,  and  9),  the  current  changes  sign. 
Since  I  do  not  allow  the  velocity  field  to  become  negative,  the  only  way  for 
the  current  to  change  sign  is  for  the  charge  density  field  to  change  sign.  But 
this  cannot  physically  happen  since  I  have  a  fluid  consisting  only  of  elec¬ 
trons.  The  current  changes  sign  because  of  the  limitations  of  the  approx¬ 
imating  functions  in  the  finite  element  expansion  (a  Gibbs  phenomenon). 
If  the  number  of  approximating  functions  is  increased,  the  frequency  of 
the  negative  current  oscillations  is  also  increased,  and  the  oscillations  de¬ 
cay  more  rapidly.  Also  note  that  the  oscillations  for  the  trapezoidal  current 
modulation  are  not  as  great  as  for  the  optimum  sinusoidal  velocity  modula¬ 
tion.  This  is  because  the  dlb/dip  entering  the  output  cavity  is  smaller  for  the 
trapezoidal  current  modulation.  One  should  not  worry  about  this  discrep¬ 
ancy,  since  the  fluid  simulation  is  only  one  step  in  a  multilayer  approach 
to  klystron  design  and  analysis.  After  achieving  a  fluid  model  optimized 
klystron,  one  should  use  the  optimized  fluid  model  as  a  starting  point  for 
the  more  detailed  particle-in-cell  (PIC)  model  of  the  electron  beam. 
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Figure  2.  Power 
convergence  for  20-kW 
beam  with  sinusoidal 
velocity  modulation. 
Equilibrium  efficiency 
is  7}  =  0.56. 


Figure  3.  Diagnostics 
for  efficiency-optimized 
buncher  cavity  with 
sinusoidal  velocity- 
modulated  20-kW 
beam.  Rf  power  is 
11.311  kW. 


Figured  Power 
convergence  for  2-MW 
beam  with  sinusoidal 
velocity  modulation. 
Equilibrium  efficiency 
is  r}  =  0.565. 
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Figure  5.  Diagnostics 
for  efficiency-optimized 
buncher  cavity  with 
sinusoidal  velocity- 
modulated  2-MW 
beam.  Rf  power  is 
1113  kW. 
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Figure  7.  Diagnostics 
for  efficiency-optimized 
buncher  cavity  with 
25  percent  duty  factor 
current-modulated 
20-kW  beam.  Rf  power 
is  17.944  kW. 
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Figure  8.  Power 
convergence  for  2-MW 
beam  with  25  percent 
duty  factor  gate 
current  modulation. 
Equilibrium  efficiency 
is  77  =  0.896. 


Figure  9.  Diagnostics 
for  efficiency-optimized 
buncher  cavity  with 
25  percent  duty  factor 
current-modulated 
2-MW  beam.  Rf  power 
is  1802  kW. 
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3.  Efficiency  Enhancement  by  Harmonic  Bunching 


3.1  Current  Modulation  for  Arbitrary  Velocity  Modulation 

To  determine  the  velocity  modulation  required  to  generate  an  efficient  cur¬ 
rent  modulation,  one  must  develop  a  way  to  calculate  the  current  modu¬ 
lation  produced  by  an  arbitrary  velocity  modulation.  One  way  is  to  use  an 
analytic  method  based  on  the  ballistic  transport  equations  (74)  and  (75): 


t 

p(z,t) 


0  U0  (to)  ’ 
Po  (to) 

1  -  zu'q/uI  ' 


and 


(74) 

(75) 


The  problem  with  this  approach  is  the  multivalued  nature  of  the  function 
that  maps  electrons  at  the  input  of  the  drift  section  to  electrons  at  the  output 
of  the  drift  section  (electrons  with  many  different  input  phases  can  have 
identical  output  phases).  Another  way  to  calculate  the  current  modulation 
that  is  simple  but  effective  is  to  divide  the  input  phases  of  the  electrons 
entering  the  drift  region  into  Nm  uniform  phase  intervals  and  the  output 
phases  into  Nout  uniform  phase  intervals  so  that 


A<pin  =  —  , 

^  jVm 

(76) 

A<pout  =  27r  , 

^  TV'out  ’ 

(77) 

<p\n  =  lApm  ,  and 

(78) 

ipf*  =  lApout . 

(79) 

In  equation  (74),  t  is  the  time  that  an  electron  requires  to  propagate  a  dis¬ 
tance  z  if  it  enters  the  system  (z  =  0)  at  time  to  arid  with  velocity  uo  (to). 
Equation  (75)  relates  the  charge  density  contribution  at  location  z  and  time 
t,  p  (z,  t),  to  the  charge  density,  po  (to),  of  the  electrons  entering  the  system 
at  time  to- 

I  can  use  equation  (74)  to  map  the  phase  grid  defined  by  equation  (78)  into  a 
grid  of  output  phases  for  a  drift  distance  of  z.  The  required  phase  mapping 
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$  (p)  is  calculated  by  renormalizing  equation  (74)  to  get  (remember  that 
ip  =  urt,  where  u>  is  the  modulation  frequency) 


In  equation  (80),  ko  is  the  wave  number  of  the  modulation  frequency  (ko  — 
uj/c),  z  is  the  drift  distance  (koz  is  dimensionless),  and  (5  (ip)  is  the  nor¬ 
malized  velocity  modulation  0  (p)  —  uo(p)/c.  Since  0  (2tt)  =  0  (0),  I  have 
$  (2n)  =  4>  (0).  Please  note  that  having  the  endpoints  of  the  interval  map 
into  themselves  does  not  imply  that  the  interior  points  of  the  input  phase 
interval  will  map  into  interior  points  of  the  output  phase  interval.  How¬ 
ever,  if  a  point  interior  to  the  input  phase  interval  maps  to  a  point  exterior 
to  the  output  phase  interval,  there  will  always  be  a  corresponding  point 
exterior  to  the  input  phase  interval  that  will  map  to  a  point  interior  to  the 
output  phase  interval.  Because  of  the  periodic  nature  of  4>  (p),  the  point  ex¬ 
terior  to  the  input  interval  will  be  related  to  the  original  point  interior  to 
the  input  interval  by  a  modulus  of  2tt.  Now  I  denote 


A$, 


=  $  (v?!n)  •  and 


(81) 

(82) 


The  relative  current  modulation  can  now  be  calculated  for  the  output  phases 
<p°ut  as  follows: 

1.  Initialize  a  relative  current  modulation  array  J°ut,  0  <  /  <  Noat  to 
zero. 

2.  Calculate  the  output  phase  array  4>/,  0  <  l  <  Nm. 

3.  Determine  if  any  of  the  output  current  density  phases  <p°ut  lie  in  the 
phase  interval  [<F/+i,  4>/]. 

4.  For  those  p°^  in  the  interval  [4>;+i.  T/],  add  Ay>In/A4>/  to  the  corre¬ 
sponding  J°ut. 

5.  If  any  part  of  the  interval  [4>/+i ,  4>/]  is  outside  the  interval  [0, 2i r),  bring 
it  inside  the  interval  [0,  27t)  by  translating  modulo  27t. 

6.  Then  test  the  translated  interval  as  in  step  3  with  any  appropriate  p°ut 
being  treated  as  in  step  4. 

7.  Repeat  steps  3  through  6  for  all  remaining  mapped  phase  intervals 

[$/+!,$/]. 


This  procedure  automatically  accounts  for  the  mapping  of  multiple  input 
phases  to  a  single  output  phase.  The  expression  A(pin/A4>/  in  step  4  results 
from  the  requirement  of  charge  conservation  in  the  system. 


Figure  10.  Current 
modulation  for  electron 
beam  with  sinusoidal 
velocity  modulation  at 
z  =  0.7,  0.8,  0.9,  and  1.0 
of  optimum  drift 
distance  zq- 


Table  2.  Parameters  for 
bunching  calculations. 


The  first  test  of  this  procedure  will  be  for  a  simple  sinusoidal  velocity  vari¬ 
ation  of  the  form 

(3  (ip)  =  (0)  (1  +  5  sin  (</>)) ,  (83) 

where  the  optimum  drift  distance  z0  satisfies  the  standard  klystron  rela¬ 
tionship  for  the  optimum  bunching  parameter  x  [3,  p  305]: 


z0ko$ 

X  =  1W 


1.84  . 


(84) 


The  results  of  the  procedure  are  shown  in  figure  10  for  the  beam  and  simu¬ 
lation  parameters  in  table  2. 

The  current  modulation  waveform  shown  in  figure  10  is  in  agreement  with 
those  shown  in  the  standard  references  [3,  p  303]. 


Variable 

Value 

E0 

35  keV 

fo 

1  GHz 

zo 

35  cm 

m 

0.37 

koZo 

7.33 

Nin 

256 

Nout 

1024 
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3.2  Optimized  Velocity  Modulation 

One  can  determine  the  optimum  velocity  modulation  waveform  for  ballis¬ 
tic  bunching  by  examining  equation  (75).  In  dimensionless  form,  the  bunch¬ 
ing  equation  is 

p{z.von{)  1 

poW")  ~  .  /?Vn)  ’  (85) 

with  <pout  given  by  equation  (80).  If  optimum  ballistic  bunching  is  defined 
as  a  repetitive  delta  function  current  modulation  at  frequency  uj,  the  re¬ 
quirement  is  that  all  input  phases  y1"  are  mapped  to  a  single  output  phase 
for  each  rf  cycle  for  a  given  z  =  zq.  This  condition  is  met  if  the  denominator 
on  the  r.h.s.  of  equation  (85)  is  zero  for  all  <pin  6  [0,  27t],  so  that  the  defining 
equation  for  optimum  velocity  modulation  is 


(87) 


To  compare  different  types  of  velocity  modulation,  I  wish  to  parameterize 
them  in  terms  of  the  average  drift  velocity  (0)  and  some  modulation  index 
6.  For  sinusoidal  modulation  (eq  (83)),  the  modulation  is  simply  character¬ 
ized  by  (0)  and  8,  where  8  =  1.84  (0)  / zqUq.  For  optimized  velocity  modula¬ 
tion,  the  situation  is  more  complicated.  I  must  integrate  equation  (87)  from 
0  to  2tt  to  compute  { 0)  and  then  solve  the  resulting  equation  for  (5  (0),  which 
is  then  back  substituted  into  equation  (76)  for  elimination.  If  one  defines  (5 
to  be 


(88) 


the  optimum  velocity  modulation  function  becomes 


(P)  1  -  e-5 


(89) 


For  the  beam  and  simulation  parameters  in  table  2,  the  result  of  the  op¬ 
timum  bunching  calculation  is  shown  in  figure  11.  The  delta  function  be¬ 
havior  of  the  optimized  velocity  modulation  as  2  ->  z0  is  clearly  shown  in 
figure  11.  The  main  problem  associated  with  the  optimum  velocity  mod¬ 
ulation  is  determining  a  way  to  produce  it.  Since  velocity  modulation  is 
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usually  imparted  to  an  electron  beam,  by  having  the  beam  transverse  a 
driven  oscillating  cavity,  one  is  limited  to  a  sinusoidal  velocity  modulation. 
One  solution  is  to  have  multiple  bunching  cavities  driven  at  frequencies  u>, 
2 u>,  3u>,  etc,  and  perform  a  Fourier  synthesis  to  approximate  the  optimum 
velocity  modulation.  Figure  12  shows  the  current  density  modulation  pro¬ 
duced  by  a  two-term  Fourier  approximation  (u  and  2 lo  modulation  terms 
present)  at  the  optimum  drift  length  for  the  exact  optimum  velocity  modu¬ 
lation.  From  figure  12,  it  is  not  obvious  why  such  a  modulation  waveform 
would  be  an  improvement  over  the  optimum  sinusoidal  waveform  shown 
in  figure  10.  One  can  argue  why  the  second  harmonic  bunching  waveform 
is  superior  by  comparing  the  two  current  modulation  waveforms  for  those 
phases  where  the  current  is  small.  This  comparison  is  made  in  figure  13, 
where  the  J  axis  has  been  expanded  to  accentuate  the  difference  between 
Jsm  and  Joptp] .  For  tp  £  [3tt/4,  5tt/4],  Jopt[2]  w  1/2  Jsin.  If  the  the  electrons  en¬ 
tering  the  klystron  output  cavity  are  classified  according  to  their  entrance 
phase,  approximately  half  the  phases  are  accelerated  and  half  are  deaccel- 
erated.  The  efficiency  of  the  output  cavity  is  critically  dependent  on  mini¬ 
mizing  the  total  kinetic  energy  (per  rf  cycle)  imparted  to  the  electrons  and 
on  maximizing  the  total  kinetic  energy  extracted  from  the  electrons.  One 
would  expect  that  during  half  an  rf  cycle,  electrons  would  be  deacceler- 
ated  and  during  the  other  half,  accelerated.  During  the  accelerated  part  of 
the  rf  cycle,  Jopt[2]  (fig- 13)  would  contain  half  as  many  electrons  as  the  Jsm 
bunching  waveform.  I  will  now  use  the  output  cavity  model  developed 
in  section  2  to  quantitatively  evaluate  the  properties  of  each  modulation 
waveform. 


Figure  11.  Current 
modulation  for  electron 
beam  with  optimized 
velocity  modulation  at 
^  =  0.7, 0.8,  0.9,  and  0.99  of 
optimum  drift  distance  zq. 
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Figure  12.  Current 
modulation  for  electron 
beam  with  second 
harmonic  approximation 
0(2]  (<£>)  to  optimized 
velocity  modulation  at 
drift  length  z0. 


Figure  13.  Comparison  of 
current  modulation  at  z0 
for  sinusoidal  velocity 
modulation  JSj„  and 
second  harmonic 
approximation  to 
optimized  velocity 
modulation  Jopt[ 2j. 


3.3  Effect  of  Harmonic  Bunching  Upon  Output  Cavity  Efficiency 

The  results  of  driving  a  simple  (single)  output  cavity  with  the  current  mod¬ 
ulation  derived  from  a  second  harmonic  approximation  to  the  optimum 
velocity  modulation  are  shown  in  figures  14  to  17.  Both  20-kW  and  2-MW 
electron  beams  are  simulated.  Figures  14  and  16  show  both  the  input  beam 
current  and  velocity  modulations  (I*,  (0)  and  /3Z  (0)  in  the  small  graph  on 
the  right  side  of  the  figures),  the  average  (over  one  rf  cycle)  rf  power  (Prf), 
and  beam  kinetic  power  (Pke),  and  the  minimum  value  of  the  velocity  field 
in  the  cavity  gap  during  the  rf  cycle,  min  (f3z).  Figures  15  and  17  show  the 
detailed  behavior  of  the  beam  current  /*>  (1)  and  beam  velocity  PZ{1)  mod¬ 
ulation  at  the  output  of  the  cavity  gap.  Additionally,  the  voltages  across 
both  the  load,  Vjoacj,  and  cavity  gap,  Vm>,  are  shown.  The  results  of  the  sim¬ 
ulations  shown  in  figures  14  to  17  are  summarized  in  table  3. 
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Note  that  since  the  second  harmonic  ballistic  bunching  is  approximate  to 
the  ideal  ballistic  velocity  modulation  shown  in  figure  11,  the  optimum 
bunching  distance  for  the  second  harmonic  approximation  may  not  be  the 
distance  implicitly  defined  by  equation  (88).  For  the  simulations  in  fig¬ 
ures  14  to  17,  the  parameter  z/zq  was  varied  in  value  near  1,  until  the  rf 
conversion  efficiency  was  maximized.  Table  3  shows  that  using  the  second 
harmonic  approximation  0[2]  (:p)  to  the  ideal  bunching  velocity  modulation 
yields  a  significant  increase  in  the  klystron  output  cavity  rf  conversion  ef¬ 
ficiency  for  both  the  20-kW  and  2-MW  electron  beams.  The  explanation  for 
this  is  evident  from  the  comparison  of  figures  15  and  3  (20-kW  beam)  or 
figures  17  and  5  (2-MW  beam).  In  comparing  the  beam  velocity  (3Z  (1)  and 
current  Ib  (1)  at  the  exit  of  the  output  cavity  for  sinusiodal  velocity  modula¬ 
tion  JSin  (figs.  3  and  5)  versus  second  harmonic  velocity  modulation  Jopt[2] 
(figs.  15  and  17),  the  major  difference  is  in  the  value  of  beam  current  outside 
of  the  main  current  bunch.  The  exit  velocity  of  the  beam  for  the  two  types 
of  modulation  is  quite  comparable  (maximum  value,  minimum  value,  and 
alignment  with  respect  to  current  bunch);  however,  the  beam  current  out¬ 
side  of  the  current  bunch  is  significantly  lower  for  the  second  harmonic 
velocity  modulation.  Since  this  is  the  phase  region  (<p)  where  the  beam  ve¬ 
locity  is  high  (for  both  modulations),  the  beam  with  second  harmonic  ve¬ 
locity  modulation  will  exit  the  output  cavity  with  significantly  less  kinetic 
power  than  the  optimum  sinusiodal  velocity  modulation.  Thus,  the  second 
harmonic  velocity  modulation  will  give  a  significantly  higher  rf . 


Figure  14.  Power 
convergence  for  20-kW 
beam  with  optimized 
second  harmonic  velocity 
modulation  and  output 
cavity  gap  located  at 
z  =  Zq.  Equilibrium 
efficiency  is  rj  =  0.713. 
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Figure  15.  Diagnostics  for 
efficiency-optimized 
buncher  cavity  with 
optimized  second 
harmonic  velocity- 
modulated  20-kW  beam 
and  output  cavity  gap 
located  at  2  =  20. 
Conserved  power  is 
(Prf)  +  (Pke)  =  20.04  kW. 


Voltage  (V) 


-3 
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Figure  16.  Power 
convergence  for  2-MW 
beam  with  optimized 
second  harmonic  velocity 
modulation  and  output 
cavity  gap  located  at 
2  —  0.92o-  Equilibrium 
efficiency  is  77  =  0.706. 


Figure  17.  Power 
convergence  for  20-kW 
beam  with  optimized 
second  harmonic  velocity 
modulation  and  output 
cavity  gap  located  at 
2  =  0.92o.  Conserved 
power  is  (Prt)  +  (Pke)  = 
2003  kW. 
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P.-(O) 


Table  3.  Summary  of 
output  cavity 
simulation  results  for 
10-/xperv  electron  beam 
with  second  harmonic 
ballistic  bunching. 


Modulation 

Beam 

power 

Cavity 

position 

Efficiency 

Efficiency 

increase* 

Power 

error* 

2nd  harmonic 

20  kW 

15  cm  (zo) 

0.713 

0.153 

+0.20% 

2nd  harmonic 

2  MW 

27  cm  (0.9zo) 

0.706 

0.141 

+0.15% 

*The  efficiency  increase  is  with  respect  to  simple  sinusoidal  velocity  modulation. 
*The  maximum  percent  value  of  ((Prf)  +  {Pke))/Pbeam  —  1-  This  value  indicates 
the  degree  to  which  fluid  simulation  does  not  conserve  energy  because  of  various 
numerical  errors. 
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4.  Conclusions 


I  have  presented  an  efficient  fluid  model  for  the  output  cavity  of  a  klystron, 
tested  the  model  for  both  realistic  and  optimized  beam  bunching,  derived 
the  optimum  velocity  modulation  for  ballistic  bunching,  derived  a  realistic 
approximation  to  optimum  velocity  modulation,  and  compared  the  oper¬ 
ation  of  a  klystron  output  cavity  for  optimized  sinusoidal  velocity  beam 
modulation  with  that  for  optimized  second  harmonic  velocity  beam  modu¬ 
lation.  I  have  shown  that  optimized  second  harmonic  velocity  beam  modu¬ 
lation  will  significantly  enhance  the  rf  conversion  efficiency  of  the  klystron 
output  cavity  at  the  cost  of  a  more  complex  beam  bunching  system  (in  the 
simplest  case,  two  driven  buncher  cavities  instead  of  one). 

In  general,  I  have  developed  an  efficient  method  of  performing  a  first  or¬ 
der  self-consistent  optimization  of  the  klystron  output  cavity.  Specifically, 
I  can  calculate  the  cavity  parameters  (for  a  given  electron  beam  and  beam 
modulation),  such  as  loaded  Q  (in  my  analysis,  I  determine  the  rf  load  and 
the  cavity  gap  gain),  that  would  maximize  the  rf  conversion  efficiency.  This 
analysis  would  be  the  starting  point  of  a  multilevel  approach  for  the  design 
of  a  klystron.  The  next  step  would  be  to  verify  the  fluid  optimization  with 
a  parametric  (cavity  modes  are  given)  PIC  code.  The  design  would  be  then 
optimized  with  the  parametric  PIC  code.  The  final  design  level  would  be 
to  simulate  the  improved  optimization  with  a  full  finite  difference  time  do¬ 
main  (FDTD)  PIC  code,  such  as  MAGIC  or  ICEPIC,  before  any  hardware  is 
fabricated. 
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